Motion of vortex sources on a plane and a sphere 
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Abstract 

The Equations of motion of vortex sources (examined earlier by Fridman 
and Polubarinova) are studied, and the problems of their being Hamilto- 
nian and integrable are discussed. A system of two vortex sources and three 
sources-sinks was examined. Their behavior was found to be regular. Quali- 
tative analysis of this system was made, and the class of Liouville integrable 
systems is considered. Particular solutions analogous to the homothetic con- 
figurations in celestial mechanics are given. 
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A point vortex model going back to G. Kirchhoff is used to describe many effects 
in classical hydrodynamics. Various problems solved with the use of this model are 
reviewed in Less known is a more general problem formulation which describes 
the motion of more complex singularities in a fluid (these singularities combine vor- 
tical properties and properties of a system of sources and sinks). This model is some- 
times of greater use for hydrometeorology. It was first obtained by A. A. Fridman, 
the famous Russian scientist who worked in hydrodynamics and cosmology, and his 
disciple P. Ya. Polubarinova (known also as Polubarinova-Kochina, April 13, 1899 - 
July 13, 1999; she was academician N. E. Kochin's wife). In Pj, the motion of two 
sources-sinks is also studied and general integrals of motion equations are examined. 
Hamiltonian equations of motion of sources-sinks were studied by A. A. Bogomolov 

(this work [3 performed independently of [2 ). The main features in the motion 
of sources and sinks (hereafter, they will be referred to as sources) are used for mod- 
eling heat convection in a fiat horizontal fluid layer, for example, periodic convection 
cells. The system of two sinks was integrated in pi; it was mentioned that the case 

> 3, as well as the case of the equations of motion of point vortices and A^ 
point masses, causes fundamental difficulties. Among the studies dealing with the 
dynamics of sources, we should also mention 01 IH], where chaotic advection of a fluid 
in a fleld of pulsating sources is examined, and a number of examples in molecular 
biology are given. 



1 Equations of motion of vortex sources 

It is well known [6J that the dynamics of point vortices in a plane (vortex lines) 
in an ideal incompressible fluid can be described by the (Hamiltonian) system of 
flrst-order equations 
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where Zj. = Xf, + iy^. are coordinates of vortices in the complex form, and F;. are 
their intensities. The velocity of fluid v = + iVy at a given point of the plane is 
determined by the expression 

1 -p 



so that the velocity near a fixed vortex z^. can be represented in the form of the 
series 

viz) = '^^ + f^a,iz-z,r. (3) 

^ n 

s=0 

Thus, the dynamics of point vortices is the dynamics of singularities in the field 
of an ideal incompressible fiuid of a special kind Q- In their work mentioned above, 
A. A. Fridman and P. Ya. Polubarinova [2] consider the dynamics of singularities of a 
more general form (these singularities persist during fiuid motion), in the case when 
the expansion of velocity near the singularity is determined by the Laurent series 

oo ^ oo 

s=l ^k) ,=0 

where a^, are complex numbers such that all do not vanish simultaneously. 
To determine the velocity of the singularity we use the following principle [2]: 

Suppose that is a ring near a singularity with the inner radius li and external 
radius I2 (see figure 1); if the fiuid bounded by the ring solidify in the 
moment considered (without redistribution of the surrounding velocity field), 
the obtained buoy will have a translational velocity vq and the angular velocity 
uoq. When the radii /i, I2 tend to zero, we obtain the respective velocities of 
the singularity. 

Using Q and taking into account the conservation laws of momentum and an- 
gular momentum, we obtain the following: 

^^0(^1,^2) = Oo, ujoih^k) = - ^2"^ ' (^) 



li + q 



and passage to the limit yields 



00, if Im 61 7^ 0, 
f = ao, uj = < (6) 
0, if Im 61 = 0. 

Thus, the translational velocity of the singularity (in the complex form) is equal 
to the absolute term in the expansion of fiow velocity into Laurent series, and the 
angular velocity does not vanish only if Im bi ^ 0, which corresponds to an ordinary 
point vortex in a fiuid with u = 00. 

Remark 1. The above reasoning shows that additional physical considerations 
(such as the existence of minimum admissible sizes within which the flow has no 



singularities and the like ), should be involved in the construction of dynamic models 
of higher- order singularities (dipoles, quadrupoles, and the like) and the evaluation 
of their angular velocity. 

The simplest singularities (which also are isotropic, that is, their angular velocity 
has no effect on their translational velocity) are the following first-order singularities 
(i.e., bi = 0,i> 2): 

1. vortex (Im hi ^ 0, Re bi = 0); 

2. sources (Im bi = 0, Re 6i 7^ 0); 

3. vortex source (Im bi 7^ 0, Re 61 7^ 0). 

The equation of vortex motion can be written as 

^| = E— . A, = J^, i=l,...,n, (7) 

f^i Zl 2-K 

where and Yi are the intensities of the source and the vortex, respectively. A 
particular case of the system ((7j) is the system of equations for motion of vortices [H] 
and sources j2l E] , which, as was found, can be represented in Hamiltonian form. 

2 Invariance of motion equations; integrals and 
Hamiltonian character 

At arbitrary complex intensities A^, equations (|7j) are invariant with respect to the 
group of motions of the plane E{^). The invariance with respect to translations 
results in the existence of two linear integrals, which in a complex form can be 
written as 

n 

Z = AkZk = const. (8) 

k=l 

In the case of ordinary vortices and sources, we have the well-known integrals 
for the vorticity center and the center of divergence: 

(for vortices) = Y1 ^kXk = const, = Y1 ^kVk = const, 

(9) 

(for sources) Qs = Yl ^k^k = const, Pg — Yl ^kVk = const. 

In the general case, Z in (jSI) determines the center of intensities. Clearly, if ^ 7^ 
0, we can choose the origin of the coordinates such that Z=0. 



Unlike the vortex dynamics, the rotational symmetry, which corresponds to a 

n 

field of symmetries of the type oi v = J2i^k^ Uk-M-), does not result, in the 

general case, in the existence of an additional integral. In this case, we have only a 
differential relation 

Y^\zki^ = Y,Kkki. (10) 

In the case of vortices ()10|) . this relationship allows us to obtain the well-known 
integral of momentum 

n 

Iv = ^ ^k{xl + yl) = const, 

k=l 

and in the case of sources, we obtain 

^ki^kVk - XkVk) = 0> ^ki^kXk + ykVk) = Y ^i^k- 

k,l 

Thus, in the case of sources, we have the non-autonomous integral 

k,l 

It is clear that, when ^^^^ ^ K^K; = 0, in the case of sources, we obtain an ordinary 
autonomous integral. In the general case of vortex sources, (jlUj) fails to yield even 
a non-autonomous integral. 

It can be easily shown that the divergence of the right-hand part of equations (|7j) 
is equal to zero; therefore, we have 

Proposition 1 The Equations of motion of vortex sources conserve the standard 
(i.e., constant- density) invariant measure. 

In addition to the above-mentioned integrals and relations, equations (|7j) in par- 
ticular cases of vortices and sources yield one more integral H and also can be 
written in the Hamiltonian form: 

ii = {xi,H}, yi = {yi,H}. 

This representation is well known for vortices; the Hamiltonian and the Poisson 
bracket are: 

i<j 

Analogous representation for sources was suggested in 3j: 

-f^ = 77- KfcK; arctan — {Xi,yA = K^^5ij. (13) 
k<l ^k — xi ■' 



The integral (jl3|) was obtained in f2]. As is the case with vortices, the integrals 
of motion Q and (jlip are non-involutive, and satisfy the following commutation 
relations: 

{Q„PJ = 5^K„ {/„QJ = -2P„ {/„Pj = 2g,. (14) 

i 

Equations (jllj) and (jl4j) . by virtue of the Liouville theorem, imply the following 
conclusion. 

Proposition 2 If in the problem of motion of three sources Kj = or ■ KjKj = 
0, the system is completely integrable in these cases. 

Indeed, if ^jKj = 0, we have three involutive integrals, for example, Q^, P,, 
and H, while if Ylij KjKj = 0, such integrals are + P^, H. 

Moreover, for the case of four sources, we obtain 

Proposition 3 Suppose that in the problem of four sources, we have Kj = 
'Yliij KjKj = 0, then the system ((7j) is completely integrable on the invariant manifold 
determined by condition = = 0. 

Indeed, in this case we have three autonomous independent integrals Q^, P^, and 
Ig, which are in involution on this manifold. Note that, for the problem of four 
vortices to be integrable, it is sufficient to require that Qs = -Ps = 0, ^ Fj = 0. 

Remark 2 In the general case of vortex sources, we can consider the complex 
function 

n = ^ AfcA; \n{zk - zi), 
which satisfies the following relations 21- 

dz^ &H -r dzj^ &H 

3 Motion of two vortex sources 

This problem is completely integrable by quadratures, and its particular cases are 
well known. 

a) Two vortices [7] 

If ri+r2 7^ 0, both vortices move along concentric circles around a center of vorticity 
with coordinates ^-^f^^f^, remaining on a common line with this center, and the 
angular rotation velocity is 



oj = where d is the distance between the vortices. 

If Fi + r2 = 0, the motion of vortices is translational in the direction perpendic- 
ular to the straight hne connecting them. And the velocity of vortices is f = 

b) Two sources El 

If Ki + K2 7^ 0, the sources move along a single straight line fixed in space and 
approach the center of intensity with coordinates ^^^f^-jj^f^; and the square of the 
distance between them varies in accordance with the following formula (P = d"^ + 
2(Ki + K2)t. 

If Ki + K2 = 0, the sources move along a fixed line so that the distance between 
them remains constant and their velocity is f = 

c) Two arbitrary vortex sources 

The equations of motion in this case are 

ll = , Z2 = . (15) 

Zi — Z2 Z\— Z2 

Suppose that Ai + A2 7^ 0; without loss of generality, we can assume that the center 
of intensity (jH)) lies in the origin of coordinates, that is 

Kizi + K2Z2 = 0. 

Expressing from here the coordinates of each particle and passing to polar coordi- 
nates = PjC*'^', i = 1,2, we obtain 

PkPk = «fc, Pk^k = -bk, k = l,2, 

where a,^ + ihf^ = {Ai + A2)~^\Ai\^, k ^ I. Solving these equations yields the law of 
motion: 

rj = Cf^ + 2a,t, = - ^ ln(Cf + 2a,t), (16) 

where (pf'^ , C^^^ are constants of integration. Thus, equations (jisp are integrable by 
quadratures. 

In this case, the trajectories of vortex sources are spirals with an infinite number 
of turns around the center of intensities. 

If Ai + A2 = 0, assuming that Z = AiZi + A2Z2, we find that 

Z2 = zi - Z/Ai, zi = — = const. 



that is, the vortex sources move with a constant velocity along parallel straight lines 
at an angle of a (tana = ■^) with respect to the segment connecting them. 

Interestingly, equations (fT^ also have an integral of the energy type, which is a 
superposition of integrals (fT^ and fT^: 

H = iI^i±lH) In - + ^ll±M. arctan ^i^. (17) 
An In x\ — X2 

It is well known that the presence of three integrals in a four- dimensional system 
(in this case, these are integrals X = Re(AiZi + A2Z2), Y = Im(AiZi + A2Z2) and H) 
implies in that such system can be represented in the Hamiltonian form [H]. If we 
take the function H in ()17j) as the Hamiltonian, the respective Poisson brackets are 



{/, g}dxi A dyi A dx2 A dy2 = cS df A dg A dX A dY, 
ae= ((Ki + K2)' + (ri + r2)')~'. 



;i8) 



Explicitly, we obtain 



{xi, Hi} = ag(K^ + r|), {xi, yj} = -seiKiKj + TiTj] 
{x„ Xj} = {yi, yj} = £e{K,Tj - K^Ti), i ^ j. 



(19) 



Remark 3 Generally speaking, the system p5|) can be represented in the Hamil- 
tonian form in three different ways, but the Poisson brackets thus obtained are not 
constant. 

Remark 4 Although the system is Hamiltonian, the merging of vortex source (at 
< 0) within finite time is not excluded. 



4 Three-sources system 
4.1 Reduction 

It is common knowledge that a three- vortex system on a plane is Liouville integrable 
and its phase space in the compact case is foliated into two-dimensional tori filled 
with quasiperiodic motions. To study that problem methods of qualitative anal- 
ysis of Hamiltonian systems (geometric interpretation, construction of bifurcation 
diagrams, and so on; see [TJ IH] , where a vast bibliography is presented) can be used. 

In the general case, the three-source problem is not Liouville integrable by 
quadratures. Its behavior, however, is regular. Indeed, this problem is determined 



by a Hamiltonian system with the Poisson bracket {xj, yj} = K,- ^Sij and the Hamil- 
tonian 

-.3 _ 

if = — V KiKAi, 9i. = arctan - — ^. (20) 
27r Xi — Xj 

Kj J 

As in the general case, the equations of motion of this system are invariant with 
respect to the E{2) group (group of motions of the plane). According to the general 
theory [S], we obtain a closed system for the invariants of the symmetry group. In 
our case (as in the vortex theory), it is convenient to take distances between the 
sources Mi = \zj — Zk\'^ , i ^ j ^ k ^ i a.s invariants. We obtain 

M, = 2( 5^ K,) + K,(M,(M-i + M-i) - M^M;^ - M,M-'). (21) 
I 

Note that such reduction into invariants of the group of symmetries in this case 
differs from the traditional Hamiltonian reduction. This is due to the fact that 
system ()20j) does not admit autonomous integrals, the Hamiltonian vector field of 
rotations is associated with the nonautonomous integral (jlip . which linearly varies 
with time. 

In the Hamiltonian form, reduction only by one degree of freedom is possible 
(based on integrals (jHl)). In this case, the distances M, should be complemented by 
one of the angles Q = 9i2] the Hamiltonian of the reduced system and the Poisson 
bracket can be written as: 



H 



M3 + Ms - Ml 
2VM2M3 



y KjKj I c) — K2K3 arccos I — I + K1K2 arccos 

yi<j J V 2VM1M3 J 

{M„ M^} = e,,kK^'A, {M3, 6} = 2{K^' + K^^^), 
{M2, 9} = K^'M^\Ms + Ml - M2), {M2, 6} = K^'M^\M3 + M2 - Mi), 

where A = Im(2;i — 2:3) (^2 ~ ^3) is the doubled oriented area of the triangle formed 
by the sources; this area can be expressed in terms of the distances by the Heron 
formula: 

= 2(MiM2 + M2M3 + M3M1) - - Ml - Ml (22) 

The equation 

• A(Mi - M2) 

® = M,M,M, ' ^^^^ 
which describes the evolution of O, completely determines the absolute dynamics 
of the sources (the remaining angles 6'23, ^13 can be found from the law of cosines). 



Thus, the Hamiltonian reduction in this case is a simple restriction on the level of 
integrals 

Let us consider equations ()21|) in more detail. Direct calculations readily show 
that system (pT|) also has the nonautonomous integral 

J2 KiK.M, -2(J2^MJ2 ) ^ = (24) 

i<j \ i<j / 

and conserves the invariant measure with the density 

p = i2A)-^ = (2M1M2 + 2M2M3 + 2M3M1 - - Ml - MlY^'^ . 

Equations (|2H) are homogenous so they also admit the one-parameter symmetry 
group Mi —>■ Mi, t —>■ At, which allows an additional reduction of the order by a 
change in the variables and the time in the following form (the so-called projection 
procedure): 

Ml = xMs, M2 = 1/M3, dt = M^xydT. 

We have 

dx 

^' = ^ = ^((Ki - ^3)y' - Ksx^ + (Ki + 2K2 + 2K,)xy - (Ki - K,)x- 

-(2Ki + 2K2 + K3)i/ + Ki), 
y' = % = ^((^2 - K3)a;2 - K3I/' + (2Ki + K2 + 2K:,)xy - (2Ki + 2K2 + K3)a;- 

-(K2-K3)i/ + K2). 

(25) 

Thus, qualitative analysis of the system considered here reduces to the examina- 
tion of the system on the plane ()25|1 . As a result, the dynamics of the system (j2n| 
is regular, although it does not have a complete set of global integrals. 



4.2 Homothetic configurations 

In vortex dynamics, an important role in the qualitative analysis of systems played 
by stationary configurations, with which stationary points of the reduced system 
were associated. Homothetic configurations (this terminology was adopted from 
celestial mechanics), i.e., configurations that remain self-similar in any time moment, 
play an analogous role in the dynamics of sources. This condition can be represented 
in the following form 

(M,Mri)-=0, z,j = 1,2,3. (26) 



Remark 5 It can be readily seen that, because of the existence of integral (I24|) . the 
system of equations (|2ip has no fixed points. 
The following simple statement is true: 

Proposition 4 The only possible homothetic configurations in a three-source 
system are: 

1. equilateral configuration, in which the sources form an equilateral triangle: 
Ml = M2 = M3; the vortices move along fixed straight lines, which cross 
in the center of intensity, and the equality 

is valid; 

2. collinear configurations, such that the vortices lie on the single fixed straight 
line Ml = z'^M^, M2 = (1 — z^M^, where z is a root of the cubic equation 

p{z) = (Ki + K2)z^ - (2Ki + K2)z^ - (K2 + 2K3)z + K2 + K3 = 0, (27) 
such that 

(Ki + K2)2(2-1)-K3 
M3 = 2 ^ ^ = const. 

z{z-l) 

PROOF. 

According to ()26|) . homothetic configurations are fixed points of the system 
The cases where x = or ?/ = correspond to singularities of the initial system, 
where one of the distances vanishes; therefore, they should be left out of consider- 
ation. The remaining fixed points can be found in the following way. Note that 
equation K3— + (Ki — K3)— = is linear with respect to y; solving this equation 

X y 

with respect to y and substituting x' = 0, we obtain, after cancellation, the equation 
that determines a homothetic configuration 

(x - 1) ((Ki + K^fx^ - (4K2 + 3Kl + 6K1K2 + 4K1K3 + 4K2K3)x2+ 

+ {3Kl + 4Kl + 4K1K3 + 4K1K2 + 6X2X3)0; - (K2 + K3)') = 0. (28) 

We find from here one more solution x = 1, which, as can be readily seen, cor- 
responds to y = 1; that is, this solution corresponds to an equilateral configuration. 
The remaining equation (PHj) . which is cubic with respect to x, with the substitution 
X = z"^ can be represented in the form 

p{z)p{-z) = 0, 



where p{z) is the polynomial (j77|) . Substituting the obtained solution for x in the 
equation for y yields y = {1 ± z)'^, that is, this solution corresponds to a coUinear 
configuration (with differently ordered sources). 

Using equation (j^Hj) . we find that in both cases, = 0, (by virtue of either 
Mi=M2 or A = 0), that is, the sources move along straight lines fixed in the 
space. The evolution rates of the squared distances can be immediately found from 
equations (jTIj) after the substitution of the obtained solutions. ■ 

4.3 Geometrical interpretation and qualitative analysis 

It is well known that complete qualitative analysis can be made for both the motion 
of a reduced system and the absolute motion p. A convenient tool for study in 
this case is a projection of the phase flow on the plane of the integral of momentum 
in the space of mutual distances Mj, i = 1,2,3. We will show that, in the case of 
three sources, such projection also exists; the corresponding plane system admits 
complete qualitative analysis and gives a complete description of both relative and 
absolute dynamics of sources. 

Let us write a nonautonomous integral of the system ()21|) in the form 

D = 2(^Ki)(^Kri)t = const, D = J2^7^Mi (29) 

and define the projection and time transformation by the formulas 

Mi , dt 
= — , dT = —, 2=1,2,3. 

The equations of motion for the projective variables are 

x'f, = 2J2^i + ^k{xk{xi^+xJ^)+x;^x,-xJ^x^)-2 (IZ^-^O (30) 

clearly, the trajectories of the system (jHUj) that correspond to the trajectories of the 
original system lie on the invariant manifold 

3 

J2^7'x, = l. (31) 

i=l 

This manifold is the plane in the three-dimensional space Xi, X2, x^. 

However, imply that the Heron formula (|22j) and the conditions Mj > 0, the do- 
main of possible motions lies within one half of the cone determined by the inequality 



2{xiX2 + X2X3 + X3X1) — xl — xl — xl > 0. 



(32) 



Thus, the boundary of the domain of possible motions on the plane (|3ip. which 
is determined by inequality (j22I), is a second-order curve of one of the three types: 



1) ellipse (if (EK,)/(KiK2K3)>0); 

2) hyperbola (if (^ K,)/(KiK2K3) < 0); 

3) parabola (if ^ Kj = 0). 

A one-to-one mapping exists between the trajectories of the system ()30|) in the 
domain considered here and the trajectories of the reduced system of three sources 
(PT|) . Thus, we have constructed a projection of the phase flow of a three-source 
system, which is completely identical to the analogous projection in the vortex 
dynamics (up to the change Kj — Fj) [Uin]- 

Let us consider in more detail the structure of the phase portrait in the first 
case (see figures 2, 3), which, by analogy with vortex dynamics will be referred to as 
"compact"; clearly, in this case all ratios Mi/Mj remain bounded for all time. The 
second ( "noncompact") case can be studied in a similar manner; however, the latter, 
third case (as well as the case ^ = 0) requires special consideration because, 
according to (j29p . an autonomous integral exists in this case. 

The structure of the phase flow of the system on the plane ()31|) is completely 
determined by its singularities, fixed points, and separatrices connecting fixed points. 
Statement 1 makes it possible to show that the system on the plane considered here 
has either four fixed points, which we denote A, Bi, B2, and (figure 2), or two 
fixed points — A and Bi (figure 3). Point A corresponds to an equilateral homothetic 
configuration and lies within the domain (ellipse); points Bi, B2, and B^ lie on the 
boundary of the domain and determine collinear homothetic configurations (since 
it is clear that for these points, A = 0). Additionally, the system has (in the 
domain considered) three singularities Ci, C2, and C3, which lie on the boundary 
(jnH) in the tangent points of planes Xi = 0; they correspond to the situation when 
the pair of sources merges. 

As in the vortex dynamics, in the compact case, we distinguish two cases, in 
which the portraits qualitatively differ: 

1°. all intensities have the same sign; 

2°. two intensities are positive (negative), while the third intensity is negative 
(positive), and Ki + K2 + K3 < 0(> 0). 



In the case of 1°, as can be seen from figure 2, coUinear configurations are saddle 
points, and tfie equilateral configuration is a node point (at Kj < 0, it is a sink, 
while at Kj > 0, it is a source), such that all the three separatrices of coUinear 
points enter into it. All the trajectories start (at Kj < 0) from singularities Ci, C2, 
and C3, do not go beyond the domain limited by the separatrices of coUinear points 
Bi, B2, S3, and merge in the fixed point A. (At Kj > 0, the direction of motion 
along these trajectories should be reversed.) 

In the case 2°, the only fixed point Bi (corresponding to the coUinear configura- 
tion) is a source (Kj < 0); the equilateral configuration now corresponds to a saddle 
point. Three separatrices of the saddle point end in the singularities Ci C2, C3, 
respectively, and the fourth separatrix, ends in the fixed point B^. A general phase 
portrait is given in figure 3. 

5 Homothetic configurations for n sources 

Of great importance in the dynamics of n point vortices are the simplest periodic 
solutions, i. e., stationary configurations such that vortices rotate around a common 
center of vorticity without changing their mutual arrangement. Analogous solu- 
tions were found in the dynamics of n point sources, with the only difference that 
the points move along fixed straight lines passing through the center of divergence 
of the sources so that their configuration remains self-similar at any time. Such 
configurations will be referred to as homothetic. 
Let us show that the following statement is true 

Proposition 5 There exists a one-to-one map between the stationary configura- 
tions of the dynamics and the homothetic configurations of the dynamics of n point 
sources. 



Indeed, for any stationary configuration in a coordinate system with the origin 
at the vorticity center, the coordinates of vortices are as follows 



where Q. is the constant angular velocity of rotation of the configuration, and the 
constant complex numbers Cfe satisfy the system of algebraic equations: 



PROOF 



k — 1., ... .,71., 




k — 1 . . .n, 



(33) 



where Ti are the intensities of vortices (A; = ^) . 

In the case of sources (A; = we similarly find the particular solutions which 
specify the homothetic configurations in the form 

Zk = Vc — 2Qt (k, k = 1 . . .n, 

where (k are complex numbers satisfying the system of equations ()33p . in which the 
change of variables Ti K/ must be made. ■ 

Note that the solution of the system of algebraic equations (jS^ is not known, 
and the problem of finding the homothetic configurations of sources also cannot 
be solved completely. A some of particular results for stationary configurations in 
vortex dynamics are summarized in p. 

Remark 6 In celestial mechanics, the homothetic configurations form a particu- 
lar case of a wider class of homographic configurations, when all points at all times 
form such configurations, which rotate around a common center of masses. There 
seems to be no such configurations, other than homothetic, in the source dynamics. 
As far as we know, the question of existence of such configurations in the case of 
vortex sources has not been considered. 



6 Vortex sources on the sphere 

Let us consider a possible generalization of the motion equations of vortex sources 
(and other singularities in a fluid flow) for the case of a sphere. To do this, we make 
a stereographic projection of the plane onto the sphere and write the equations of 
vortex motion on the sphere pQ in complex form. 

In the case of the stereographic projection from the southern pole onto the plane 
passing through the northern pole (see figure 4), the rectangular coordinates of a 
point on the sphere can be written as 

, = ( , = a , ji i-M.e+,f) 

where A = (2i?)~^ is the curvature. 

Using the complex coordinates of vortices (k = + iVk, k = 1 ... n on the plane 
considered, the motion equations of the vortices can be represented in the form 

^ is " IS M?) j ■ '''' 



Remark 7 The factor before the second summand in (|35|) in the hook J7|/ appears 
to he erroneous. 

The velocity of the fluid determined by the vortices in any point of the sphere is 
expressed analogous to (jH^j) : 

' ^ is " IS n M?) j ■ '''' 

It can be easily shown that the divergence of this velocity field is zero and its vorticity 
is constant at any point of the sphere C 7^ Cfc and equals 



k=l 



Y.^k. (37) 



By analogy with the planar case, let us replace the intensities iVk by complex 
numbers = + iVk] we obtain the equations describing the motion of sin- 
gularities. These singularities can be interpreted as vortex sources on the sphere. 

-^(5A-(S"')^)- - 

Unlike the planar case, the velocity field determined by the singularities (this 
field can be calculated by formula (jHBjl with replacement iVi — > A^), in addition to 
the vortex sources in points Cfc, has a non-zero vorticity (jHTj) and a divergence equal 
to 

1 " 
1=1 

at any point of the sphere C 7^ Cfc- 

Clearly, such corrections in the motion equations of a fluid on a sphere stem from 
the compactness of the sphere, and if a point source or a vortex appears in some 
point, at the same time, a point sink or sinks uniformly distributed over the surface 
of the sphere must appear somewhere. Equations have not been studied yet. 
It is still unclear whether they have at least one first integral and whether they are 
Hamiltonian. 
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Figure 1: 




a) Ki = K2 = K3=l h) Ki=l, K2=0.5, K3 = 0.125 

Figure 2: Trajectories of the system ()3n|l on the plane (j3T|) in the case 1°. 




Figure 3: Trajectories of the system (jHUj) on the plane (jHTj) in the case 2° (Ki = 1, 
K2 = 0.5, K3 = 0.125). 





N q 








/e/2 J 


s 





Figure 4: 



Figure 1. 

Figure 2. Trajectories of the system (jHDI) on the plane (j!?T|l in the case 1°. 
Figure 3. Trajectories of the system (jHUI) on the plane in the case 2° (Ki = 1, 
K2 = 0.5, K3 = 0.125). 
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